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ABSTRACT 

Context. A recent modeling study of brightness ratios for CO rotational transitions in gas typical of the diffuse ISM by Liszt found the 
role of H collisions to be more important than previously assumed. This conclusion was based on recent quantum scattering calcula- 
tions using the so-called WKS potential energy surface (PES) which reported a large cross section for the important — > 1 rotational 
transition. This result is in contradiction to one obtained using the earlier "BBH" PES for which the cross section is quite small and 
which is consistent with an expected homonuclear-like propensity for even AJ transitions. 

Aims. To revisit this contradication with new scattering calculations using two new ab initio PESs that focus on the important long- 
range behavior and to explore the validity of the apparent departure from the expected even AJ propensity in H-CO rotational excita- 
tion obtained with the WKS PES. 

Methods. Close-coupling (CC) rigid-rotor calculations for CO(v = 0, J = 0) excitation by H are performed on four different PESs. 
Two of the PESs are obtained in this work using state-of-the-art quantum chemistry techniques at the CCSD(T) and MRCI levels of 
theory. 

Results. Cross sections for the J = — > 1, as well as other odd AJ, transitions are significantly suppressed compared to even A J 
transitions in thermal energy CC calculations using the CCSD(T) and MRCI surfaces. This is consistent with the expected even AJ 
propensity and in contrast to CC calculations using the WKS PES which predict a dominating — » 1 transition. 
Conclusions. Inelastic collision cross section calculations are sensitive to fine details in the anisotropic components of the PES and 
its long-range behavior. The current results obtained with new surfaces for H-CO scattering suggest that the original astrophysical 
assumption that excitation of CO by H 2 dominates the kinetics of CO in diffuse ISM gas is likely to remain valid. 
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1. Introduction 

Carbon monoxide is found in a variety of astrophysical environ- 
ments usually being the most abundant molecule with the ex- 
ception of H2. It is typically observed in absorption in the UV 
and near-IR and in emission in the far-IR to millimeter wave- 
lengths. For example, the Infrared Space Observatory (ISO) 
detected a large number CO rovibrational lines toward Orion 
Peak 1 and 2 (Gonzalex-Alfonzo et al. 20021) while pure rota- 
tional lines were observed in the planetary nebulae NGC 7027 
dCernicharo et al. 1997b . 

In most environments, the molecular level populations are 
not in equilibrium, so that predicting spectral line intensities de- 
pends crucially on the magnitude of collisional excitation rate 
coefficients due to the dominant species H2, He, and H. Adopting 
semiempirical potentials, Green & Thaddeus (1976) performed 
one of the earliest quantal scattering calculations for the CO- 



H2, CO-He, and CO-H systems at astrophysically relevant tem- 
peratures. Their results suggested that the rate coefficients for 
atomic hydrogen collisions were significantly less than those due 
to H2, so that CO-H collisions were usually neglected in model- 
ing studies of CO spectra. 

The H-CO complex has been extensively studied by the 
chemical physics community for the past three decades includ- 
ing early ab initio PES calculations by Botschwinna (1974). A 
later ab initio PES by Bowman, Bittman, and Harding (1986, 
BBH) was constructed and used in subsequent scattering cal- 
culations (Lee & Bowm an 19871) . The H+CO system has a bar- 
rier to the formation of the HCO complex at a H-CO internu- 
clear separation of 4.2 ao with a calculated height of 0.125-0.17 
eV, but measured to be 0.087(17) eV (-700 cm- 1 or 1000 K) 



( |Wang et al. 1973[ ). Lee & Bowman (1987) found that for col- 
lision energies below the barrier, and not coincident with tun- 
neling resonances, rotational excitation cross sections displayed 
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an even AJ propensity, i.e. homonuclear-like behavior at long- 
range. 

Another H-CO surface was computed at the MRCI level of 
theory by Keller et al. (1996), the so-called WKS PES. This PES 
was found to reproduce H-CO spectroscopic measurements bet- 
ter than that of the BBH surface. It was also noted by Green et al. 
(1996) that the WKS and BBH surfaces gave different rotational 
excitation cross sections with the difference being most signif- 
icant for the — > 1 transition. This discrepancy was ascribed 
to the differences in the treatment of the long-range portion of 
the potentials. Balakrishnan, Yan, & Dalgarno (2002) later con- 
firmed this difference in the rotational excitation cross sections 
and therefore utilized the WKS surface for large-scale calcula- 
tions of rotational excitation rate coefficients for astrophysical 
modeling applications. 

Recently, Liszt (2006) adopted the new CO-H rotational 
excitation rate coefficients of Balakrishnan et al. (2002) in an 
investigation of the non-equilibrium rotational populations of 
CO in environments typical of the diffuse interstellar medium 
(ISM). He concluded that the increase in the rotational excitation 
rate coefficients for CO-H substantially enhanced CO rotational 
brightness and excitation temperatures and that CO-H collisions 
should not be neglected in future models. 

However, what has been previously overlooked is that the 
new large — > 1 rate coefficient results in a departure from the 
homonuclear-like even AJ propensity which has been argued 
on reasonable physical grounds to be valid for the H-CO col- 
lision system JLee & Bowm an 19871 1. As the scattering results 
have been confirmed in a number of theoretical investigations, 
a possible explanation is a subtle deficiency in the long-range 
portion of the WKS surface. In this letter, we explore this is- 
sue by computing two new and completely independent H-CO 
rigid-rotor surfaces and performing close-coupling calculations 
to obtain the rotational excitation cross sections. We briefly de- 
scribe the potential surface computations and scattering calcu- 
lation methods in the next two sections. The results and their 
astrophysical implications are presented in Section 4 followed 
by a conclusion in Section 5. 



2. Potential Energy Surface Calculations 

One of the two new ab initio surfaces was calculated using 
the coupled cluster method with single and double excitations 
and a perturbative treatment of triple excitations [CCSD(T)] 
employing the frozen core approximation (Purvis & Bartlett 
1982, Raghavachari et al. 1989). For the open-shell calcula- 
tions, restricted open-shell Hartree-Fock orbitals were used in 
the CCSD(T) calculations, but the spin-restrictions were relaxed 
in the solution of the coupled cluster equations [R/UCCSD(T)] 
(Knowles et al. 1994). We adopted the doubly augmented cor- 
relation consistent d-aug-cc-pVnZ (n = T, Q, and 5) basis set of 
Woon & Dunning (1994). The full counterpoise correction was 
applied to approximately account for basis set superposition er- 
ror (Boys & Bernardi 1970). At each configuration, calculations 
were performed with all three basis sets for the dimer (HCO) 
and monomers (CO and H) and these energies were extrapo- 
lated to the complete basis set (CBS) limit. These extrapolated 
total energies were then used to compute the CBS limit coun- 
terpoise corrected interaction energies. The interaction potential 
was computed in Jacobi coordinates with the CO bond length r 
fixed at its experimental equilibrium value (2.1322 ao) (Huber 
& Herzberg 1979). The resulting rigid-rotor PES that is a func- 



tion of R (distance between H and the CO center of mass) and 6 
(angle between the vectors r and R) was fit to the expression 

A' 

V{R,ff) = Y J P^0)V x {R) (1) 

where the are normalized Legendre polynomials with N = 
12. The Va(R) were calculated using 12th-orderGauss-Legendre 
quadrature. The two-dimensional grid on which the ab initio 
calculations were carried out corresponded to 23 values of R 
between 3-20 ao and the 12 values of 6 were obtained from 
the 12 Gauss-Legendre quadrature points. The large range in 
R was chosen in order to obtain good coverage of the long- 
range interaction, which is important in the low-energy scatter- 
ing, as discussed below. The calculations were carried out with 
the MOLPRO suite of ab initio programs dAmos et al. 2 002). 

The second PES was computed using the com- 
plete active space self-consistent field [CASSCF] 
(IWerner & Knowles 19821 IKnowles & Werner 19851 ) 

and internally contracted multireference configuration 
interaction [MRCI] ( IKnowles & Werner 19881 [19921 
IWerner & Knowles 19891 1 method with aug-cc-pVQZ basis 
set for H, C, and O by Woon & Dunning ( 1994) as implemented 
in the MOLPRO program. Calculations were performed for 
1-3 1 A' states in the C s point group for all non-linear approaches 
in Jacobi coordinates, and for a CO bond length of r — 2.20 ao. 
The rigid rotor PES was obtained for the following grid of 
points: R = 1.6-5.0 (0.2), 5.0-5.8 (0.4), and 6.0-15.0 a (1.0 a ) 
and 0=1-179° (15°), where the numbers in parentheses indicate 
the step-sizes for the corresponding range and coordinate. 
The final energies incorporate the Davidson correction ( 119741 1. 
The computed potential energy values were combined with an 
accurate ab initio long-range potential that includes Ce, Oj, 
and Cs coefficients for the H-CO interaction. These coefficients 
have been obtained using Pade approximants with Cauchy 
moments evaluated by linear response coupled cluster (LR- 
CCSD) theory (Hattig et al. 1997) with the d-aug-cc-pwCV5Z 
basis set (Peterson & Dunning 2002) and correlating all 14 
electrons of CO in the CCSD calculation. The long range part 
has been smoothly joined to the ab-initio calculated points and 
a full analytical representation of the rigid rotor potential with 
accurate long-range interaction was constructed and used in the 
scattering calculations. 



3. Collisional Excitation Calculations 

In the present study, the quantum close-coupling (CC) method 
(Arthurs & Dalgarno 1960, Lester 1971, Secrest 1979) was used 
for the calculations of rotational excitation due to H atom im- 
pact. CO was treated in the rigid-rotor approximation being held 
fixed at its equilibrium bond-length. The CC scattering calcula- 
tions were performed by expanding the angular dependence of 
the relevant potential in Legendre polynomials as given in Eq. 
(1). Calculations were performed with the non-reactive molec- 
ular scattering progam MOLSCAT ( IHutson & Green 19941 to 
generate integral state-to-state cross sections. The CC equations 
were integrated using the modified log-derivative Airy propa- 
gator of Alexander and Manolopoulos d 19871 1 with a variable 
step-size. A sufficient number of total angular momentum partial 
waves and basis set terms have been included to secure conver- 
gence of the cross sections. 
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4. Results and Discussion 

Figure la compares the first four terms (A = - 3) of the 
Legendre expansion for the BBH and WKS surfaces as well 
as the new PESs computed in the current work, referred to as 
CCSD(T) and MRCI. The agreement between all surfaces ap- 
pears to be reasonable at the scale of the figure. However, focus- 
ing on the long-range portion, clear differences are apparent in 
the A = 1 component as shown in Fig. lb. The present PESs are 
believed to be more reliable than previous surfaces because we 
have specifically focused on the long-range behavior through ex- 
plicit calculations at much larger internuclear distances than in 
the earlier work. Further, the importance of the long-range be- 
havior and its influence on low-energy collisions is elevated due 
to the barrier at R = 4.2 ao as described below. 

Calculations of state-to-state cross sections from initial rota- 
tional state J = for the excitation of CO by H on the MRCI, 
CCSD(T), BBH, and WKS surfaces were performed for colli- 
sion energies at 400 and 800 crrT 1 as given in Figs. 2 and 3, 
respectively. The cross sections using the BBH and WKS sur- 
faces are in good agreement with those obtained by Green et al. 
(1995) and Balakrishnan et al. (2002), respectively. While the 
results presented here adopted the rigid-rotor approximation, we 
performed tests with BBH and WKS and found that the result- 
ing cross sections using the full surfaces, i.e., with vibrational 
motion, are practically identical. The largest difference was less 
than ~10% for the — » 1 transition with WKS. Our rigid -rotor 
cross sections in Figs. 2 and 3 further confirm the previously 
noted result that the — > 1 transition using WKS is approx- 
imately an order of magnitude larger than that obtained with 
BBH. 

Comparing now the rigid-rotor cross sections obtained with 
the new CCSD(T) and MRCI surfaces, we see that i) the two 
sets of cross sections are in very good agreement, ii) they re- 
produce the even AJ propensity obtained with BBH, and iii) the 
— > 1 transition is small with a value similar to that obtained 
with BBH, in contrast to the WKS result. These observations can 
be interpreted by examining the A = 1 terms displayed in Fig. 
lb since the dominant contributions to the state-to-state cross 
sections originate from potential coupling matrix elements with 
A = \AJ\. The A - I components from the CCSD(T) and MRCI 
calculations are seen to be very similar, both displaying a peak 
near 5 ao and comparable long-range behavior. Conversely, for 
the WKS surface the A = 1 component peaks at a shorter inter- 
nuclear distance, closer to the barrier location, and has a different 
long-range behavior. Indeed, the long-range tail of the WKS po- 
tential appears to be decaying slower than the other potentials. 
The BBH A = 1 component is smaller than that of the other sur- 
faces which may explain its slightly smaller — » 1 cross section. 

To further test the apparent discrepancy related to the A = 1 
component of the Legendre expansion, we examined scattering 
results for excitation of the J = 1 and 2 rotational states of 
CO(v = 0). The cross sections for the AJ-1 excitation transi- 
tions with the WKS surface were found to also be an order of 
magnitude larger than those obtained with the other surfaces. 

Based on all of our scattering calculations in conjuction with 
the computation of the new surfaces, it appears that the WKS 
surface has a different anisotropic behavior, especially at long- 
range. This anistropy, evident in the A — 1 component, while 
of small magnitude, drives the AJ = 1 transitions and results in 
the significant discrepancy in the important J = — > 1 cross 
section. 

Our findings have important implications for modeling of 
CO in non-equilibrium situations, e.g. diffuse interstellar gas, 
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Fig. 2. Excitation cross sections due to H impact for CO(v = 
0, J = 0) to CO(v' = 0, J') for a center of mass collision energy 
of 400 cm -1 . Results obtained with different PESs are indicated 
by the symbols. 




Fig. 3. Same as Fig. 2, but for a center of mass collision energy 
of 800 cm- 1 . 



where atomic hydrogen is present. The conclusions of Liszt 
(2006), based on the adoption of the rotational excitation rate 
coefficients obtained with WKS, that H collisions can signifi- 
cantly affect the CO rotational populations, and therefore model- 
ing of the resulting absorption/emission lines, is no longer valid. 
In particular, the — > 1 cross section due to H impact is found 
to be much smaller than that due to H2 collisions in agreement 
with the earlier results of Chu & Dalgarno (1975) and Green & 
Thaddeus (1976). Further, our results illustrate the sensitive de- 
pendence that low-energy scattering can have on fine details of 
the potential energy surface. 

5. Conclusions 

The current scattering calculations using the new CCSD(T) and 
MRCI surfaces have corroborated the expected even AJ propen- 
sity for the H-CO system at low collision energies. These find- 
ings indicate that the rotational excitation rate coefficients ob- 
tained in Balakrishnan et al. (2002) using the WKS surface are 
inaccurate bringing into question any astrophysical conclusions 
deduced from modeling with these rate coefficients. Instead of 
being large, the J = — * 1 cross section is small as originally 
found by Chu & Dalgarno (1975). The supposition of Green & 
Thaddeus (1976) that excitation due to molecular hydrogen and 
He should dominate the non-equilibrium kinetics of CO remains 
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valid. Calculations are currently underway to compute the full 
long-range HCO surface at the CCSD(T) level which will allow 
for new large-scale scattering calculations for an extensive range 
of rotational and vibrational transitions. 
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